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Abstract — The problem of minimizing the total power con- 
sumption while satisfying different quality-of-service (QoS) re- 
quirements in a two-hop multiple-input multiple-output network 
with a single non-regenerative relay is considered. As shown 
by Y. Rong in |1|, the optimal processing matrices for both 
linear and non-linear transceiver architectures lead to the di- 
agonalization of the source-relay-destination channel so that the 
power minimization problem reduces to properly allocating the 
available power over the established links. Unfortunately, finding 
the solution of this problem is numerically difficult as it is not in 
a convex form. To overcome this difficulty, existing solutions rely 
on the computation of upper- and lower-bounds that are hard to 
obtain or require the relaxation of the QoS constraints. In this 
work, a novel approach is devised for both linear and non-linear 
transceiver architectures, which allows to closely approximate the 
solutions of the non-convex power allocation problems with those 
of convex ones easy to compute in closed-form by means of multi- 
step procedures of reduced complexity. Computer simulations are 
used to assess the performance of the proposed approach and to 
make comparisons with alternatives. 

Index Terms — MIMO, non-regenerative relay, power alloca- 
tion, majorization theory, transceiver optimization, quality-of- 
service requirements, non-convex optimization, power consump- 
tion, closed-form solution, decision-feedback equalizer. 



I. Introduction 

THE demand for high-speed and high-quality multime- 
dia services in wireless communication systems has 
increased significantly over the last years. This has led to 
a strong interest in multiple-input multiple-output (MIMO) 
technologies as they represent a promising solution to improve 
the system reliability and coverage by means of space-time 
coding technique |2j and/or to increase the spectral efficiency 
through spatial multiplexing fS) and I?]. The main impair- 
ment of MIMO systems is represented by the multi-stream 
interference (MSI) arising from the simultaneous transmis- 
sion of parallel data streams over the same frequency band. 
However, if channel knowledge is available at the transmitter 
and receiver side appropriate linear or non-linear transceiver 
architectures can be used to mitigate MSI and fully exploit 
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the potential benefits of MIMO technologies iQ. The design of 
transceiver architectures for MSI mitigation in MIMO systems 
has received great attention in the last years and several 
excellent works can be found in the open literature. Good 
surveys of the results obtained in this area can be found in 
161 - ID in which the authors develop general frameworks for 
the optimization of linear and non-linear MIMO architectures. 
This is achieved by resorting to the notions of additive and 
multiplicative majorization theory (see [|9] for a complete 
reference on majorization). 

In the case of a long transmitter-receiver distance, single 
or multiple repeater (relay) nodes may be necessary to pass 
information from the transmitter (source) to the receveir 
(destination). The re-transmission schemes employed at the 
relays may operate according to several different protocols. 
Among them, the decode-and-forward protocol makes use of 
a regenerative relay to first decode the received signal and 
then re-encode and forward the original information to the 
destination node. On the other hand, the amplify-and-forward 
protocol adopts a non-regenerative relay in which the received 
signal is first linearly processed and then re-transmitted toward 
the destination. Clearly, amplify-and-forward is more suited to 
practical implementation than decode-and-forward since de- 
coding and re-encoding multiple data streams involves higher 
computational complexity and larger processing latency than 
simply amplifying and forwarding them. Among the different 
amplify-and-forward MIMO relay systems, the simple two- 
hop single-relay MIMO model has gained a lot of interest 
over the last years as it provides a reasonable tradeoff between 
potential benefits and practical implementation issues. For 
the above reasons, a two-hop single-relay MIMO network is 
considered in this work. 

The optimization of MIMO non-regenerative relay networks 
has received much attention recently. A survey of the results 
obtained in this area can be found in |ld| and briefly summa- 
rized in the following. One of the first attempt in this direction 
can be found in fTTl and ||T2| in which the relay amplifying 
matrix is designed so as to maximize the link capacity between 
source and destination. It turns out that optimizing the relay 
matrix largely improves the system capacity with respect to 
alternative solutions based on heuristic arguments. Although 
the capacity is one of the most important information-theoretic 
measure, there are many other ways of characterizing the 
reliability of transmission schemes. For this purpose, several 
different solutions based on linear processing in the form 
of decorrelating or minimum mean-square-error schemes as 
well as on non-linear layered architectures have been recently 
derived and investigated according to different criteria. A 
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unifying framework for the design of linear transceivers in 
the presence of a single relay is presented in llT3l and later 
extended to the multiple relay case in lfT4l . The main results 
of ifTSll and lfT4l are achieved by application of additive 
majorization and consist in proving that the joint optimization 
of source, relays and destination matrices under fixed power 
constraints diagonalizes the MIMO relay channel as long as 
the objective function is Schur-concave or Schur-convex. In 
the latter case, the diagonalizing structure is optimal provided 
that the transmitted data symbols are properly rotated before 
channel diagonalization. The results in fT3\ and |fT4l are of 
great interest since many different optimization criteria driving 
the design of communication systems arise in connection 
with Schur-concave or Schur-convex functions [ITSI . All the 
above results have been recently extended in llT6l to the 
case in which a decision feedback equalizer (DFE) is used 
at the destination node. This is achieved by means of the 
multiplicative majorization theory and the equal-diagonal QR 
decomposition tool illustrated in ifTTl . In particular, it is found 
that when the objective function is multiplicatively Schur- 
convex the optimal design leads to the uniform decomposition 
of the MIMO relay channel into an arbitrary number of 
identical parallel single-input single-output relay subchannels. 
This leads to a transmission scheme characterized by a much 
lower bit-error-rate (BER) than a system employing linear 
processing at the destination. On the other hand, if the objec- 
tive function is multiplicatively Schur-concave the optimum 
non-linear transceiver architecture reduces to the linear one 
discussed in |[T3l and | fT4l for which the channel diagonalizing 
structure is optimal. 

All the aforementioned works are focused on the minimiza- 
tion or maximization of a global objective function subject to 
fixed power constraints at source and relay nodes. This may 
prevent their use in those multimedia applications supporting 
several types of services each characterized by a different 
reliability constraint. A common approach to overcome this 
problem is to deal with the minimization of the power con- 
sumption while meeting the quality-of-service (QoS) require- 
ments for each data stream (see for example llTSi - ||201 and 
references therein). A solution in this direction is proposed 
in II2TI in which source and relay matrices are designed so 
as to ensure a specific signal-to-noise ratio (SNR) on each 
subchannel (see also 1221 for single-antenna relay networks). 
In in, the author makes use of majorization theory and 
propose a unifying framework for the design of linear and non- 
linear transceiver architectures that minimize the total power 
consumption either in single-hop or multi-hop MIMO relay 
systems. As in lfT3l - lfT6l . it turns out that the optimal solution 
leads to the diagonalization of the source-relay-destination 
channel. However, the resulting power allocation problem is 
only upper- and lower-bounded using a successive geometric 
programming approach and a dual decomposition technique, 
respectively. Unfortunately, the computation complexity of 
both solutions is relatively high so as to make them unsuited 
for practical implementation. An alternative solution with 
reduced complexity is proposed in |!23l where the authors rely 
on a convex relaxation of the QoS constraints. This results into 
two convex suboptimal problems that provide different upper 



and lower bounds to the optimal solution. Unfortunately, such 
bounds do not meet each other over the entire region of interest 
but only for high values of SNR. This means that they cannot 
be used to exactly characterize the solution of the original 
problem. Furthermore, it is worth observing that relaxing the 
QoS constraints may result into a suboptimal solution which 
does not necessarily belong to the feasible set of the original 
problem. This is the case of the lower bound illustrated in |l23| 
in which a rescaling operation is required to meet the original 
QoS constraints, thereby leading to an increase of the power 
consumption. 

In this work, we return to the problem of designing optimal 
linear and non-linear transceiver architectures in a two-hop 
single-relay MIMO network with QoS requirements and we 
make use of the theoretical analysis presented in yj to focus 
only on the resulting non-convex power allocation problems. 
The latter are tackled following a different approach that 
provides a framework in which the power allocation for linear 
and non-linear schemes can be considered in a unified and 
concise way, and allows us to extend the results given in fTTl. 
In particular, the non-convex power allocation problem is first 
reduced to an equivalent form and then approximated with a 
convex one, whose solution is within the same feasible set 
of the original problem and can be computed in closed-form 
through a multi-step procedure that requires no more than K 
steps (where K is the number of data streams). Numerical 
results are used to highlight the effectiveness of the proposed 
approach. Interestingly, it turns out that for both linear and 
non-linear systems the approximated solutions are very close 
to those of the original problems. 

The remainder of this paper is organized as foUowsQ. Next 
section describes the two-hop system model and introduces 
the power minimization problem. In Section III, a linear 
transceiver architecture is considered and the proposed ap- 
proximated solution is described together with the convexity 
analysis of the original power allocation problem. In Section 
IV, the results are extended to a non-linear architecture in 
which a DFE is used at the receiver Simulation results are 
discussed in Section V while some conclusions are drawn in 
Section VI. 

II. System description and problem formulation 

We consider a flat-fading MIMO network in which the 
information data are carried from source to destination with 

'The following notation is used thi'oughout the paper. Boldface upper and 
lower-case letters denote matrices and vectors, respectively, while lower-case 
letters denote scalars. We use A = diag{a,i ; n = 1, 2, . . . , K} to indicate 
a K X K diagonal matrix with entries a„ while A^'^ and A^/^ denote the 
inverse and square-root of a matrix A. We use Ik to denote the identity matrix 
of order K and rank(A) to indicate the rank of a matiix A while [ J^. ^ is the 
(fc, £)th entiy of the enclosed matrix. In addition, we use E {■} for expectation, 
I for the Euclidean nonn of the enclosed vector and the superscript and 
for transposition and Hermitian transposition, respectively. The notation 
< xA-y > stands for xy = 0, x > and y > 0. If the elements of x 
and y are arranged in increasing order, we use x ^+(™) y to say that x is 
weakly additively majorized by y. 

^Although specific for a flat-fading channel, the model adopted throughout 
the paper can easily be extended to frequency selective environments using 
orthogonal frequency-division multiplexing (OFDM) as a transmission tech- 
nique. 
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the aid of a single non-regenerative relay. The information 
first flows from source to relay and then from relay to 
destination. The direct link between source and destination 
is not considered as it is assumed to undergo relatively large 
attenuation compared to the link via the relay. 

The fcth symbol is denoted by and is taken from an 
L— ary quadrature amplitude modulation (QAM) constellation 
with average power normalized to unity for convenience. We 
denote by K the total number of transmitted symbols and 
assume that source and destination are equipped with N 
antennas while the relay has M antennafl 

The source vector s = [si, S2, ■ ■ ■ , sk]'^ is first linearly 
processed by an iV x if matrix U and then transmitted over the 
source-relay link in the first time-slot. At the relay, the received 
signal is processed by an M x M matrix F and forwarded 
to the destination node in the second time-slot. The received 
signal at the destination takes the form lfT3l 

r = HUs + n 

where H = H2FH1 is the equivalent channel matrix while 
Hi and H2 are the source-relay and relay-destination chan- 
nel matrices, respectively. In addition, n is a zero-mean 
Gaussian vector with covariance matrix pRn, where Rn = 
H2FF^H|'^+Iivf and p > accounting for the noise variance 
over both linko Henceforth, we denote by 



and 



Ho 



^H2 ^ H2 



the singular value decompositions (SVDs) of Hi and H2. 
Without loss of generality, in all subsequent derivations we 
assume that the entries of the diagonal matrices and 
Ah^ are arranged in non-increasing order This amounts to 
saying that Ah,,„ > A^i.n+i and Xn^.n > ^H2,n+i for 
71 = 1,2,..., A' - 1. 

The vector r is finally processed by the destination node 
for data recovery. For this purpose, we consider either a linear 
receiver or a non-linear detector equipped with a DFE. In both 
cases, we deal with the joint design of the processing matrices 
so as to minimize the total power consumption given by 

Pt = tr {UU^} + tr {F (HiUU^Hi^ + pl^) F^} 

while satisfying different QoS requirements in terms of the 
MSEs. As mentioned before, the above problem has recently 
been discussed by Y. Rong in Ul. In the next, the main results 
of HI are first briefly revised and then the major contributions 
of this work are described. The linear case is considered first. 

III. Linear Transceiver Design 

When a linear K x N receiver G is employed, the vector 
y at the input of the decision device takes the form 

y = GHUs + On. 

^The results can be easily extended to a more general case in which different 
number of antennas are available at source and destination. 

''The extension to the case in which the noise contribution over each link 
has a different variance is straightforward. 



In these circumstances, the constrained power minimization 
problem can be mathematically formalized as follows JT]: 



mill Pt 

U.F.G 



s.t. 



[E]„,„ < Vn 



for n= 1,2,...,^' (1) 



where E denotes the MSB matrix E = E{(y — s) (y — s)^} 
while the quantities < < 1 are design parameters that 
are assumed to be in non-decreasing order, i.e., ?7„ < rjn+i- 

As proven in 11], under the assumption that K < 
min(rank(Hi), rank(H2)) the optimal G^"^ in ([T]i is equal 
to the Wiener filter while the optimal U'°^ and F'°^ have the 
following form 



H 



and F^°>=Vh.AJ^ 



'Hi 



(2) 



where V^f^, V/f^ and CIhi correspond to the K columns 
of , and CIhi associated to the K largest singular 
values of the corresponding channel matrix while Q is a 
suitable K x K unitary matrix such that 



Vn 



for n = 1, 2, 



The matrices Ajj ~ diag{A[j^j, Ajjj' ■ 
diag{A^°', , A^ 



Ao) x(o) 
^_F,1' -F,2' 

entries given by 



..,K. 
X'i^]^} and Af 



(3) 



^fk} ^^^^ diagonal structures with 



A 



(°) 

U,n 



-Ai°^ and A 



On 



(4) 



where An^ 
problem: 



and B 



(o) 



' " >^H2,n Ai°'> + 1 

are the solutions of the following 



K 



mm 

{A„ >0},{S„ >o} 



n=l 



Ar. 



B„ 



A_f/i,n Xh2, 



(5) 



s.t. ELi An < ELi Vn for i = 1, 2, . . . , if 
0<A„<1 n=l,2,...,K 
A„<A„+i n= l,2,...,ii:- 1. 

Here, A„ denotes the nth eigenvalue of E and is given by fl] 

An+Bn + l 



1 



(6) 



From dill, it is seen that U*-"^ and F^"-* are obtained as 
matched filters along the singular vectors of Hi and H2, 
respectively. In addition, the overall channel matrix given by 
n = G(°)H2F(°)HiU(°) reduces to U = QA^Q" where 
is diagonal with entries |1J 

\(o) \„ \(o) \„ 



1 



As discussed in HI, the above result reveals that the optimal 
structure of the relay communication system is diagonal up to 
a unitary matrix Q satisfying (O. If K is a power of two, such 
a matrix can be chosen equal to the discrete Fourier transform 
matrix or to a Walsh-Hadamard matrix llTSl . Otherwise, it can 
be determined through the iterative procedure described in 
|2l. 

The only problem left is to solve (|5]l or, equivalently, to 
properly allocating the available power on the established 
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links. This is a challenging task since the quantities A„ in (|6]l 
are not jointly convex in An and i?„, thereby resulting into 
a combinatorial minimization problem for which no practical 
algorithm is available. A way out to this problem consists in 
solving (|6]l alternately with respect to A„ and _B„ keeping the 
other fixed. This leads to an iterative optimization procedure 
that if properly initialized monotonically converges to a local 
optimum of (|6]l since the conditional updates of An and _B„ 
may either decrease or maintain (but not increase) the objective 
function. Although conceptually simple, the above approach 
requires to iteratively solve multiple convex problems and does 
not guarantee the convergence to the optimum. An alternative 
approach is discussed in [1] in which the optimal solution 
is upper- and lower-bounded using a geometric programming 
approach and a dual decomposition technique, respectively. 
Unfortunately, the computation complexity of both solutions 
is relatively high so as to make them unsuited for practical 
implementation. To overcome the above problems, we follow 
an alternative approach in which the optimization is first 
carried out over An and Bn for a fixed A„ and then over all 
possible A„ within the feasible set of ©. In Appendix A, using 
standard calculus techniques it is shown that the quantities An 
and Bn that minimize the function p{An/ XHi,n + Bn/ \H2,n) 
in (|5]l for a fixed A„ have the form 



of Ai°^ and 



and 



An 



Bn — 



1 — A„ / Xhi^u %/! — A„ 



An y ^Hi.n A„ 

Using (Q and (O into p{An/ Xhi.u + Bn/XH2,n) yields 

Pn (An 



(7) 



(8) 



P 



\/XHi,nXH2 

where 7„ is given by 



^J_K^^VTZ^^ (9) 



A„ 



A. 



XHi,n + XH2,n 



(10) 



,n 

It is worth observing that 7„ is a positive function of XHi,n 
and XH2,n taking values in the interval [2,oo). The minimum 
7„ = 2 is achieved for Xn^.n = XH2,n- 

The optimization over all possible A„ satisfying the con- 
straints in (|5]l leads to the following equivalent problem: 

K 

min P(A) = ^ (A„) (11) 



where A = [Ai, A2, . . . , A^]^ and C is the set of admissible 
A defined as 

£ = |a e R^' : ^+('^' A, < A„ < 1 and A„ < A„+i} 

with T] = [771, 772, . . . , 77a']^ being the vector collecting the 
QoS requirements. 

To proceed further, we denote by A*^"^ = 
{x[°\xi°\...,x'-°'^} the solution of (HB and call P(A(°') 
the corresponding power consumption. Once computed, the 
nth element of A'°^ is used in ^ and dHJ for the computation 



which are then employed in (|4]i for 
determining A[^^„ and A^^„. 

Unfortunately, finding A'"-* is a challenging task since the 
optimization problem in (fTTT i is not in a convex form ||251 . 
While the feasible set C is convex, the function P„(A„) is 
convex for < A„ < a„ and concave for a„ < A„ < 1 
where, as shown in Appendix B, a„ is given by (see Fig. 1 
for a graphical illustration) 

_ ^ (12) 

3 
4 



with ^„ defined as follows 

^""4" ^16(^2 _ 4)+ ^4(^2 _ 4)2- (13) 

The two above expressions hold for 7„ 7^ 2. For 7„ = 2, a„ 
takes its minimum value given by (see Appendix B) 



y 

This result has an interesting theoretical relevance as it allows 
to prove the following lemma. 
Lemma 1: If 




K 

E 

n=l 



< 



(14) 



then the optimization problem in ( fTTT i is convex. 

Proof: If (O holds true, then J2n=i ^» < 8/9 for all 
A G £ from which (bearing in mind that A„ > 0) it follows 
that any admissible A„ must be smaller than or equal to 8/9. 
Using this fact and recalling that P„(A„) is convex in (0, 8/9], 
we have that the objective function P(A) is convex VA G C 
as it is the sum of convex functions. ■ 

Lemma 1 establishes a sufficient condition for the optimiza- 
tion problem (fTTT i to be convex. Clearly, such a condition is 
not always met as it depends on the number of data streams 
and on the specific QoS requirements. This means that solving 
(fTTT i is in general hard and prompts us to search for alternative 
methods. As a major contribution of this work, in the next 
section (fTTT i is approximated with a convex problem whose 
solution is close to the optimal one and can be evaluated in 
closed-form through an exact procedure requiring a maximum 
number of K — 1 steps. 

A. Hyperbola-based approximation 

We start observing that -Pn(A„) in (|9]l has approximately a 
1/A„— shape over the interval (0,a„] in which it is convex. 
This holds true especially when A,i approaches zero. Then, the 
idea is to approximate P„ (A„) with the following hyperbola 
(see Fig. 1): 

Hyp„ (A„) 



Wn 
An 



where u)„ and z„ can be obtained as follows 



arg rnin sup 

•"•'^ A„e(o,i] 



Pn {Xn) 



An 



(15) 
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Fig. 1. Graphical illustration of Pn{\n), iiyp„ (^n) and Ln{\n) for < 
A„ < 1. 



It can be shown that Wn and z„ are given b}0 

^ : (7n + 2) 



and 



p 27„ + 3 



(16) 



(17) 



,71 

Replacing P„(A„) with Hyp„ (A„) in (fTTT i leads to the follow- 
ing convex optimization problem: 



mm 

AG £ 



where we have omitted the irrelevant terms {z„}. 

The above problem is clearly in a convex form as the 
objective function and the feasible set are both convex. To 
solve it, we first observe that the the ordering constraint 
An < A„+i in C is always satisfied by the solution of ( fTsT l. 
To see how this comes about, recall that we have assumed 
Ahi,„ > \hi,7i+i and Xh^^u > Ans^n+i- Then, from ^ 
it follows that Wn < u'n+i- Now, denote by A* the solution 
of ( fTsT i and assume that it is such that A* > A*_|_]^ for some 
n e {1,2,..., A' — 1}. Then, the vector A' obtained from 
A* by simply exchanging A* and X^+i would still satisfy 
the convex constraints of ( fTSI l but it would result in a lower 
objective function since w„ < Wn+i- This means that A* must 
be such that A* < A*^^, for n = 1, 2, . . . , A' - 1. 

The numerical evaluation of the solution of (fTsT l cannot 
be performed by means of classical bisection or hypothesis 
testing methods but it requires the development of specific 
algorithms exploiting the particular structure of the problem. A 

^The proof has been omitted for space limitations. It will be provided upon 
request. An intuitive explanation at least for w„ relies on the observation 
that P„(A„) in (9) goes to infinity as p {-fn + 2) / (A„ ^XHi,n>^H2,n) for 
An — > 0. This means that when A,i approaches zero the supremum in l |15t 
is bounded only if the coefficients w„ are in the form given by U6) . 



Algorithm 1 Multi-step procedure for solving (IT¥t . 
1) Set i = 0, Sq = Q and compute 



6 J = ?7„ 



forj = l,2,...,if. 

2) While i < K: Compute 



^K^i = min(l, Ak_j) 



where 

Xk-i 

Set 



max 



=0,1 K-i-1 y^"-' rrjy- 



Sj for j = 0,l,...,A:-i-2 

Sk-i - Xk-i for j = K - i-l 



and i = i + I. 



good sample in this direction is represented by the multi-level 
waterfilling algorithm illustrated in lfT9l (see also ll26l ) which 
provides the solution of (fTSI l in no more than K{K + l)/2 
iterations. In this work, we propose an alternative algorithm 
that allows the computation of A* in closed-form with a 
maximum number of K steps. 

Proposition 1: The solution to (fTsT l can be computed 
through the iterative procedure illustrated in Algorithm 1. 

Proof: See Appendix C. ■ 
As seen from Algorithm 1, differently from lfT9l a closed- 
form solution depending on the system parameters is provided 
for each step of the iterative procedure. It is worth observing 
that a similar algorithm has been recently proposed in 11271 
for the optimization of linear two-hop MIMO networks with 
multiple relays. Although similar, the proposed one has been 
derived following a different line of reasoning which is not 
specific for the optimization problem at hand but it can be 
used to easily accommodate other optimization problems with 
similar structures. As we will see in the next section, it can be 
applied to a MIMO relay network in which a DFE is employed 
at the destination to compute an approximation of the optimal 
power allocation. 

Once A* is computed through Algorithm 1, it is then 



used for approximating Ajjjj and A]?„ in (|4]i using the same 
procedure illustrated before for A'^°^ More precisely. A* is 
first employed in (|7]i and ^ to obtain A* and S*, which are 
then used to replace An '' and Bn '' in (01). This yields 



U,7l \ 1 \* 



1 ~ -^n , /Agi,,i y/l - A* 



and 



^F,7l 



^/g^(l-A*) + VT^ 



while the required power results equal to P{\*). 
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It is worth observing that replacing A*-"^ with A"^ inevitably 
increases the power consumption. Then, we have that 

P(A*) > P(A(°)) 

with P{\^°'^) being the global-minimum power allocation for 
satisfying the QoS requirements. The numerical evaluation 
of the performance loss P{\*) — P(A'°-') would require 
knowledge of P(A'^°^), which can be computed only solving 
(fTTT i. A possible way out to this problem is to make use of 
a lower bound of P(A^°''), say P. More precisely, observing 
that P < P(A(°') and P(A(°') < P(A*) yields 

P(A*)-P(A(°') <P(A*)-P. 

The above result indicates that the performance loss is upper- 
bounded by the difference P(A*) — P. Clearly, if such a 
difference is relatively small, the power penalty incurred 
in approximating A'"-* with A* can be neglected and A* 
well approximates the solution of the original problem in 
(fTTT i. Motivated by the above fact, we now proceed with the 
computation of a lower bound for P(A^°^). To this end, we 
start observing that for any A„ e (0, 1] P„ (A„) is not smaller 
than the following convex function; 



where 



Pri ('^Ti ) 7 



< A„ < f3n 



PnWn) + PLWn) {Xn - f3n) , /3„ < A„ < 1 



7„ + 4^) (19) 



denotes the first derivative of P„(A„) evaluated at A„ — /3„. 
The latter represents the abscissa of the intersection point 
between P„ (A„) and the Une P,' (A„ — 1) passing through 
(1,0) and tangent to P„(A„) (see Fig. 1). This amounts to 
saying that 

PniPn) + P'nWn) [K - Pn) = P'nWn) (A„ - 1) 

or, equivalently, that P!^{(3n){Pn — 1) = PniPn) from which 
using (|9]l we obtain 



Pn{Pn)Pn 



In Vl - A. + 2 = 0. (20) 



Substituting (dUl into ^ yields 3^„ - 2 = 7^(1 - /3„)^- It 
turns out that the solution of this equation in the interval (0, 1] 
is unique and given by 

Pn^l- , (21) 



[en + 1)^ 



with 



3/72 -2 + 7„ ^7,2 -4 



(22) 



Now, consider the following optimization problem: 

K 

min L(A) = ^L„(A„) (23) 



obtained by replacing P„(A,i) with Ln{Xn) in ( fTTT ). The above 
problem is in a convex form since both the objective function 



and the feasible set are convex. Moreover, denoting by A = 
{All A21 • • ■ 7 Aa:} its solution we may write 

L(A) < P(A(°)) 

where we have used the fact that A'^"'' belong to the same fea- 
sible set of ( I23] ). In addition, as a consequence of the following 
inequality L(A) = ^fL, L„ (A„) < Pn (An) = P(A) 

we have that 

p(a(°)) < p(a(°)). 

Collecting the above facts togheter yields 

L(A) < i(A(°)) < P(A(°') 

from which it is seen that a lower bound for P(A*^°^) is given 
by P = L{\). Numerical results shown later demonstrate that 
the difference P(A*) — P is negligible. As discussed before, 
this makes A* a very good approximation of A'"-*. 

Following a simple line of reasoning, it can easily be shown 
that P represents the tightest lower bound of P(A'^°') that 
can be obtained by replacing P(A) in ( fTTT i with a convex 
lower approximation. To see how this comes about, take a 
look at Fig. 1 and observe that L„(A,i) represents the convex 
function not greater than Pn{Xn) for which the difference 
Pn{Xn) — Ln{Xn) IS minimum for any A„ G (0, 1]. This makes 
Ln{Xn) the best convex lower approximation to Pn(A„) in 
the interval (0, 1]. Accordingly, L{X) is the best convex lower 
approximation to P(A) in (0, 1]^^. Since P = L(A), the above 
statement follows easily. 

IV. Non-linear Transceiver Design 

When a non-linear receiver with DFE is employed, the 
vector y at the input of the decision device (under the 
assumption of correct previous decisions) can be written as 

y = (GHU-B)s + Gn 

where B is a strictly upper triangular matrix B of order K. 
The power minimization problem is formulated as: 

min Pt s.t. [E]„ „ < for n = 1, 2, . . . , A' (24) 

where Pt is given by (HJ while E denotes the MSB matrix. 

As proven in U, the optimal G^°-' in (|24] | is the Wiener 
filter while B(°) takes the form B(°) = DL^-I^. The matrix 
L is lower triangular and such that 

LL" = u(°)"h(°)"rL°)"h(°)u(°) + piK 

with H(°) = H2F(°)Hi and rL°' = H2F(°)f(°)''h|^ + Im 
while D is diagonal and designed so as to scale to unity the 
entries [DL^]„^„ for n = 1, 2, . . . , A'. The optimal U^") and 
F(°) take the form 

V'-"^ ^Vh.J^u^S" and F^"' = V^^ A^/^fig, 

where S is a suitable K x K unitary matrix such that 

[L];'„ = forn=l,2,...,i^. (25) 
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The entries of the diagonal matrices Au and Aj? are still given 
by © with the quantities and si"'' now obtained as: 



500 



mm 

{A„ > 0}.{B„ > 0} 



K 

n=l 



B„ 



Hi,n 



A 



H2,n 



(26) 



s-t. Ili=i>^n<Ili=iVn for J- = 1,2, 
0<A„<1 n = l,2,...,A' 
A„ < A,i+i n 1,2, . . . , A' - 1. 

As for the linear case, the optimal processing matrices lead to 
a channel-diagonalizing structure provided that the symbols 
are properly rotated at the source and destination nodes by 
a unitary matrix S chosen such that dZSl l is satisfied. This is 
achieved through the iterative algorithm illustrated in |l28l. 

Following the same procedure as in the linear case, we first 
minimize over An and Bn for a fixed A„ and then over all 
possible A„ within the feasible set of ( |26] |. The first step yields 

K 

min P(A) = Vp„(A„) 



(27) 



n=l 

s-t. riLi^n <nLi^" for j = i,2, 

< A„ < 1 n = 1,2, . . . , A' 
A„<A„+i n = l,2,...,A'-l. 



Observe that the feasible set of the above problem is not 
convex since the inequality constraints 

n' 1 ^ n' 1 J = 1, 2, ... , K (28) 

do not form a convex set for AT > 2 125)1 . Letting 

6'n = In A„ and k„ = In ?7„ 
the problem can be readily reformulated as follows 

K 

min Q(0) = yQ„(0„) (29) 
e e A/ — ' 

n— 1 

where = [9i, 6*2, ... , ^k]"^ and Qn (On) takes the form 







1 



(30) 



,? 

while A/^ is defined as 

A" = 1^ e : K -<+^"') 6>,e'„ < Oand6'„ < 6'„+i} 

with K = K2, . . . , Kif]^. 

Although A/^ is convex, the problem (|29] l is not in a convex 
form due to the properties of function Qn{On)- The latter 
is shown in Appendix D to be convex for 0„ < liKpn and 
concave for In 0„ < 6'„ < (see Fig. 2) where 0„ is given by 

0n = 1 , ^ = (31) 



with 



8 ^ 



Xn = 3 + ^ V277^i ' 104+ V(277,^ - 104)2 -16+ 



104 - ^(277^ - 104)2 



16. 
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Fig. 2. Graphical illustration of Qn{6n). Exp^{9n) and for —1.5 < 

e„ < 0. 



The minimum value of (pn is achieved for 7„ = 2 and given 
by (see Appendix D) 



2(V2-1 



(32) 



It is worth observing that the above result cannot be used as 
done in the previous section to derive a sufficient condition on 
the QoS requirements under which the convexity of (|29^ can 
be established. 

As for the linear case, in the next ( |29] l is replaced with a 
convex problem easy to solve whose solution is shown to be 
close to the optimal one. 

A. Exponential-based approximation 

Using the same methodology illustrated in the previous 
section, Q„(6'„) is replaced by the following exponential 
function: 



Wn 



Exp„(6l„) 

where u)„ and z„ are obtained as follows 

(w„, z„) = argmin sup Q,i (6'„ 

"''^ e,.e(-oo,o] 

Using standard analysis not shown for space limitations, it is 
found that w„ and z„ are still given by (fT6] l and (fTTl i. 

Then, the non-convex power allocation problem in (|29] | is 
approximated with the following convex one: 



K 

71 — 1 



(33) 



whose solution 6* = [9\, 62, . ■ ■ , 0'^]'^ is easily proven to be 
such that 9n < S^+i- Interestingly, it turns out that 6* can be 
computed using an iterative procedure derived following the 
same line of reasoning adopted to obtain Algorithm 1. 
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Algorithm 2 Multi-step procedure for solving d3Jt . 

1) Set z = 0, (5o = and compute 

3 3 
^3 ^^Kn^^\n ?]n 
n—1 n—1 

forj = l,2,...,if. 

2) While i < K: Compute 

e*ji_^ = min(l,^K_,) 

where 



Ox-i = lnu)A'-i+ max 

e=0,l,...,K-i-l 



^K-i - Si - J2f=i+i In'' 



Set 

63 = 
and i = i + 1 



K~i~l 
for j = 0,l,...,X-i-2 



' 1 6K-^-0l._.^ iOT]=K-l~l 



Proposition 2: The solution to ( |33] | can be obtained through 
the iterative procedure illustrated in Algorithm 2. 

Proof: See Appendix E. ■ 

Once 6* is computed through Algorithm 2, it is then used to 
approximate a[^^„ and A^|^ exactly in the same way illustrated 
for the linear case. This produces 



and 



\* — 



1 - e}>r. 



-(1 



A 



where we have used the fact that A* = e^". 

In order to validate the quality of 6*, a lower bound is now 
computed. Paralleling the steps of the previous section, the 
lower bound is obtained as Q = S{6_) where 6^ is the solution 
of the following problem: 

K 



min 5(0) = ^5„(0„) 



n=l 



(34) 

0n < In Ipn 

(ln?/'„) + Q^j(ln'0„)(6'„ - InV'rO , Ini/'n < 07i < 

2-iJr 



As shown in Fig. 2, Sn {On) is given by 
where 

Qn (InV'ri) = ] , ' , I In + 



1pn\/^Hi.n^H; 



2.n \ VI - V'ti 

evaluated at 0„ = In'f/'n 



is the first derivative of Qn{ 
which represents the abscissa of the intersection point between 
Qn (On) and the line passing through the origin and tangent 
to Qn{&n)- In Appendix D, it shown that ■4>n is such that 

In i^n ( , 2 - ■(/)„ 



In 



VI - i'n 



= 0. (35) 




Fig. 3. Total power consumption wlien equal QoS constraints are given with 
TV = 3 and p = 1 or 0.01. 

TABLE I 

Total power consumption when equal QoS constraints are 

GIVEN WITH N = 3 AND p = 1. 



Algorithm 


r? = 0.9 


r? = 0.5 


Tj = 0.1 


r? = 0.05 


r? = 0.01 


L-HA 


1.001 


14.211 


28.5907 


32.020 


39.316 


L-GP 


1.109 


14.254 


28.612 


32.121 


39.321 


L-LB 


0.990 


14.113 


28.5907 


32.020 


39.316 


NL-EA 


0.356 


11.291 


23.351 


27.248 


34.960 


NL-GP 


0.361 


11.314 


23.359 


27.250 


34.981 


NL-LB 


0.325 


11.252 


23.330 


27.247 


34.959 



Unfortunately, the solution of the above equation cannot 
be computed in closed-form but it can only be evaluated 
numerically. 

As for the linear case, numerical results shown later demon- 
strate that the difference between Q{6*) and Q is negUgible. 
Moreover, Q still represents the best convex lower approxi- 
mation of the original problem in (|29] | . 

V. Numerical results 

Numerical results are now given to highlight the effective- 
ness of the proposed solutions. Comparisons are made with the 
successive geometric programming (GP) approach illustrated 
gin m. The CVX convex optimization toolbox for MATLAB 
is used to solve the optimization problem. The number of 
antennas employed at the source, relay and destination nodes 
is equal, i.e., N = M. The number of transmitted symbols is 
fixed to K = A^. The entries of Hi and H2 are modeled as 
independent complex circularly symmetric Gaussian random 
variable with zero mean and variance 1/iV. All numerical re- 
sults are obtained averaging over 10^ independent realizations 
of the channels. 

Fig. 3 illustrates the total power consumption as a function 
of the QoS constraints when iV = 3 and the noise variance 
is p = 1 or 0.01. The same QoS constraint is imposed for 
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Fig. 4. Total power consumption when equal QoS constraints are given with 
TV = 4 and p = 1. 
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Fig. 5. Total power consumption when unequal QoS constraints are given 
with N = i and p = 1. 



each data stream, i.e., rjn = ij for n = 1,2, . . . , K. The curves 
labelled with L-HA and NL-EA refer to a system in which 
the power is allocated using Algorithm 1 and Algorithm 
2, respectively. On the other hand, L-GP and NL-GP refer 
to a system in which the successive GP approach of is 
employed in conjunction with a linear or a nonlinear receiver, 
respectively. The results of Fig. 3 indicate that the total power 
consumption of the proposed solutions is substantially the 
same as that obtained with the schemes discussed in |[T|. 
However, this is achieved with much lower complexity as 
L-HA and NL-EA are obtained in closed form by means 
of Algorithm 1 and 2. Moreover, both algorithms attain the 
corresponding lower bounds denoted by L-LB and NL-LB and 



evaluated solving ( |23] | and ( |34] |. To ease comparisons, some 
of the results of Fig. 3 are also listed in Table I] As already 
shown in JT], it turns out that the non-linear architecture 
provides the best performance for all the investigated values of 
?/. As expected, increasing the noise variance from 0.01 to 1 
enhances the required power of all the investigated solutions 
of approximately 20 dB. Similar conclusions can be drawn 
from the results of Fig. 4 in which = 4 and p = 1. 

The results of Fig. 5 are obtained in the same operating 
conditions of Fig. 4 except that now rji = rj/A, 772 =773 = ?7/2 
and ri4 = rj and p = 1. Compared to the results of Fig. 4, the 
total power consumption increases due to the more stringent 
requirements over some established links. 

VI. Conclusions 

We have discussed the power allocation problem in a two- 
hop MIMO relay network in which the goal is the min- 
imization of the total power consumption while satisfying 
different QoS requirements given in terms of the MSEs. 
The original non-convex power allocation problem has been 
approximated for linear and non-linear architectures with a 
convex one that can be solved exactly through a multi-step 
procedure of reduced complexity. Comparisons with existing 
alternatives requiring much higher computational burden show 
that the same total power consumption is required. It is worth 
observing that the extension of the proposed approach to the 
case in which multi-hops are used to carry the information 
from the source to the destination is not simple and it is 
currently under investigation. 

Appendix A 
Keeping A„ fixed, from ^ we have that 



1 



An + A„^„ — 1 

while p{A,JXHi.n + Bn/XH^^n) reduces to 

Br, 



(36) 



P 



Xhi,ji 

An , 1 



P 



An 



- 1 



_XHi,n XH2,n \Xn + A,iA„ — 1 

Taking the derivative with respect to An and equating it to 
zero yields 



. _ 1 — A„ / Xhi ,n a/1 — A„ 

An — ^ ± 



^^2 



from which using ( |36] | yields 



„ 1 — A„ / A/f2,n \/l ~ -^ri 

i>ri = ^ ± 



An y XE_\,n Xn 

We now observe that the solution 



An — 



1 — Xn / A//i,„ \J\ — A„ 



An 



XH2,n Xn 



and 



Bn = 



1 — Xn / XH2,n VI — Xn 



Xn 



X 



Xn 



(37) 



(38) 
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is not admissible as it violates one of the constraints An > 
and Bn > in Q. To see how this comes about, note that 
An in ( [37] i is larger than zero only if XHi,n < A^f2.„(l — A„). 
On the other hand, from (l3Ft we have that Bn > only if 
Aifa.n < — A„). Collecting these two facts together, 

it is seen that the constraints An > and i?„ > are both 
satisfied when AHi,n/(l - A„) < XH2,n < Ah,,„(1 - A„). 
Since < A„ < 1 it follows that such inequalities cannot be 
ensured by any value of \Hi,n and \H2,n- Then, the optimal 
An and i?„ are given by (Q and (|8]i, respectively. 

Appendix B 

The convexity of the objective function Pn{\i) in ( fTTT i on 
V = {A„; < A„ < 1} is studied. For notational convenience, 
the index n is omitted. 

We start taking the derivative of P (A) in (|9]l with respect 
to A and obtain 

P' (A) = , -(1+ 4=4") (39) 



\/l~A 



from which it easily follows that P'{\) < for all A £ {VXl}. 
This means that P (A) is a decreasing function of A. Taking 
the derivative of P' (A) yields 



P"(A) = 



/(A) - 7,9(A) 



(40) 



2AVAHiAff, (l-A)VT^ 

where we have defined /(A) = 3A^ — 12A + 8 and 3(7) = 
-4(1 - A) Vl - A. Setting /(A) = 7g(A) leads to the fol- 
lowing quartic equation 

9A'* - 8 (9 - 27^) A^ + (4 - 7^) {ASX^ - 48A + 16) = 0. 

After lengthy computations (not shown for space limitations) 
it turns out that there exists an unique point A in solving the 
above equation. For 7 7^ 2, such a point A = a is given by (fT2] i 
in the text. On the other hand, setting 7 = 2 into the quartic 
equation produces A^(9 — 8A) = from which it follows that 
a = 8/9. Once a has been computed through (fT2] i and ( fT3] l, 
we have that /(A) — 7(?(A) > for any < A < a. This 
amounts to saying that P" (A) > or, equivalently, that P (A) 
is convex over the convex set < A < a. Vice versa, for any 
a < A < 1 we have /(A) - 7g(A) < so that P" (A) < 
and the function P (A) is concave. 

A close inspection of (fTZb and (fTST l reveals that a depends 
exclusively on 7. We are now interested in computing its 
minimum value a,„in as 7 varies. Although amin could be in 
principle computed using the closed-form expression given in 
(fTZt . a more simple line of reasoning is followed henceforth. 
Observe first that /(A) in ( l40l i is a convex and decreasing 
function while g{X) is a negative, concave and increasing 
function in V. Moreover, recall that 7 is a positive function of 
and A^^ taking values in the interval [2, 00). Collecting 
these facts together, it follows that the point a such that 
/(a) = 7.9(a) moves towards the upper limit of as 7 
increases. This means that its maximum value is given by 
ttmax = 1 and it is achieved when 7 goes to infinity. On the 
other hand, amin is achieved when 7 takes its minimum value, 
i.e., 7 = 2. As shown before, this yields amin = 8/9. 



Appendix C 

In the next, we highlight the major steps leading to the 
solution of (fTSl l in the form given by Algorithm 1. For this 
purpose, we adopt the following approach. We first compute 
\k, which represents the largest solution of the problem 
obtained from ( fTsT i after removing the constraints A„ < 1 
for n = 1,2,..., A'. Interestingly, Xk can be obtained in 
closed-form and it can be efficiently used to compute A^. 
All the above results are then used to derive a simple iter- 
ative procedure providing all the remaining solutions A* for 
n = 1, 2, . . . , AT — 1 in no more than AT — 1 iterations. 

The first step requires to find the largest solution of the 
following problem: 



mm 

{0<A„} 



S-t- Ei=l An < 



1,2,.. 



where 



(41) 



(42) 



To this end, we denote by "P^*) for i = 0, 1, 2, . . . , AT - 1 
the optimization problem obtained from (l4ll by removing the 
inequality constraints from the (i + l)th to the (K — l)th, i.e.. 



Wn 

mm > , — 

{0<A„} ^"=1 A„ 

s-t. Ei=l An < Sj 
J2n=l An < 



J =l,2,...,i 



from which it follows that ^'(^^i) is equivalent to i4T\ while 
V^^"^ is obtained from (l4TT l by removing all the constraints 
except the last. Since the above problem satisfies the Slater's 
condition, the nth solution of V^^K say Xn \ is found from the 
necessary and sufficient Karush-Kuhn-Tucker (KKT) optimal- 
ity conditions. This yields 



w„ 



7— n ^' 



(0 



for n = 1,2, 



and 




for n = i + 1, j + 2, . 



where the Lagrange multipliers <;j are chosen to satisfy the 
following constraints: 

0<.f±(^.-ELAi'')>0 (45) 

for j = 1,2, . . . ,i and j — K. From ( l44l i. it is found that 



,A' 



(43) 



(44) 



a(* 



A 



K 



for n = i + 1, i + 2, . . . , A'. (46) 



Moreover, it is easily seen that ^j^^ > from which using (l45T l 
it follows that the last constraint is always satisfied with strict 
equality, i.e., An^ = (dis- 

similarly, the solutions of p'^+i) take the form 



w„ 



y-J+l (i+1) 



for n = 1,2, . . 



,i+l 
(47) 
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and 



ioTn^i + 2,i + 3,...,K (48) 



from which it follows that the solutions A}^ for any i < K — I 
can be easily computed once A^'' is given. The latter is easily 
obtained from ( fSSl l and reads 



where the Lagrange multipliers ^j'^'^^ are such that 

o<^-^)±(^,-5:^^^A(:-))>o 



Ai°) 



'JK 



EA' , 
1=1 



for n^l,2,...,K 



(49) 



for j ~ 1,2, . . . ,i + 1 and j = K. 

To proceed further, we focus on 7^(*+i) and consider the two 
conditions <;f|+^^ > and <;f|+^^ = 0, separately. If "J^+t^^ > 0, 
from (|49]l it follows that the (i + l)th constraint of T'f^+i) is 
satisfied with equality, i.e., 

AW =<5,+i. (50) 

^ — ^n— 1 

Also, since ^n^^'^ — we may write 

^ — ^n—t-\-2 ^ — ^71—1 ^ — ^n— 1 

(51) 

Recalling yields A^i+^^ = A^+^' ^ for n = i + 
2, . . . , iiT from which we may write 



where we have defined do = 0. Using the above result and 
applying repeatedly dSSl l, it turns out that the solution of 
interest given by Xk 
form as follows 



can be determined in closed- 



/WK- 



5t+i 



:+2 



In addition, it is easily recognized that if J2n=i ^" 
then 

E /.''>'^^+i 

^ — ^n— 1 



(52) 



so that 



Substituting (|46] | into the above inequality produces 



\ (j) ^ / ^ ^'+1 



(53) 



ri=i+2 



Putting (|52]) and ([531), it is found that when > then 



A 



a: 




On the other hand, if ^^^^^^^ = from (|49] l we have 

^•^+l 



■"^^ — ^n— 1 



Accordingly, we may write 

E'' > - -^^'+1 



and 



Xk >V^^ —■ 



(54) 



In this case, it is easily seen from (|43] l - ( |45] l and Wf\ - ( |49] | 
that T^^*) and T'^^+i) have the same solution. In particular, this 
means that A^^ ~ ^k^^^- 

Collecting the above results togheter, it turns out that 



max I A]^ , y/wK 



5k - 5i+i 



(55) 



n=i+2 V , 



Afc- = \/wk max — — 



Once \k has been computed through 
determined by applying the following lemma. 
Lemma 2: The solution A^ is always such that 



(56) 



A^- can be 



Xk 



miii(l, Ak)- 



(57) 



Proof: Using the KKT optimality conditions of (TSl) it is 
seen that 



A„ — 



fern = 1,2,..., if 



(58) 



where the Lagrange multipliers {Cj } and {7,1} must be chosen 
to satisfy the following constraints: 



< 



0±U,-E„^i^« >0 j = l,2,...,X (59) 



and 



0<7„±(A;-1)>0 71 = 1,2,..., 



(60) 



In order to prove (|57t , the two cases A^^ < 1 and A^^ = 1 are 
considered separately. Assume first AJ^ < 1. Then, from ( I6OI 1 
we have 7^ = 0. In addition, since the solution of ( fTSb is such 
that A;^ < X*^j^i for n = 1, 2, . . . , K~l, it follows that A;'^ < 1 



and 7„ = Vti. Accordingly, ( 158b - (1591 1 become formally 
equivalent to the KKT conditions of problem ( |4T]) (see ( l43b 
- ( l45l ) with i = K — V) meaning that the two problems have 
the same solution. In particular, A^ = \k ~ min(l,A/f). 
Assume now X*j^ = 1. In this case, it is easily recognized that 
Xk > 1- Indeed, if by absurd Aa' < 1 then the solution of 
( I4T] ) would coincide with the solution of ( fTSl ). and A^ would 
be equal to Xk (less than 1), which contradicts the hypothesis 
X*K = 1- Then, even in this case X^ = 1 = min(l,AA') and 
the result in dSTl i easily follows. ■ 
To see how all the above results can be used to find the 
solution of ( fTSI ). observe that when A^ has been computed 
through ( l56l l and (l57l i. the remaining solutions A* for n ~ 
1, 2, . . . , A' — 1 can be found by solving the following problem 



mm > , 

{0<A,.<1} ^"=1 A„ 

E,i=i Xn <Sk - X*i^. 



s.t. 



-I Xn < 5j 



j = 1,2,...,K -1 
(61) 



Recalling that 6j = Eri=i write 



Sk - Xk 



^ * n — 1 
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Assume now A^, = 1. Observing that < ?7„ < 1 we have 



Sk - A^- < Sk^i. 

The same result can be obtained when At 
([56] l. This means that ( |6TI ) reduces to: 

„A'-1 Wn 

{o<A„<i} '^"-i A„ 



K 



Xk < I using 



s.t. 



ELiA«<<5^ j-l,...,A--l 



(62) 



where 



- A}, 



for j l,2,...,X-2 

for j = K — 1 . 



We now proceed computing AJ^^ j^. For this purpose, we 
follow the same procedure used for the computation of 
We first compute \k-i, which is the largest solution of the 
optimization problem obtained from (l62t after removing the 
constraints A„ < 1 for n ~ 1, 2, . . . , A' — 1. Paralleling the 
same steps leading to (|56] | yields 



\k-1 = ^WK-l 



from which A5f_]^ is obtained as 

A^r-i = ™n (1, Ak-i). 

Applying repeatedly the same steps for any n < K—2 leads to 
the iterative procedure illustrated in Algorithm 1 in the text. 

Appendix D 

The convexity of the objective function Q„(0„) in ( l30l l on 
Q = {6'„;0ri < 0} is studied. To simplify the notation, we 
omit the index n. Taking the first derivative of Q (6) produces 



Q'{0) 



\/ Ah2 

from which it follows that 



7 + 



(63) 



Q" (0) = 



P 



7 + 



6e« 



2(1 



.3/2 



Setting Q" (9) = and letting 



yields 



-2j{l — (f))"^/^ or, equivalently, taking the square of both sides 

(3 - 7^) </)3+4 (11 - 37^) (1)^+4: (4 - 7^) (1 - 30) ^ 0. 

It turns out that the above equation has a unique solution in 
the interval (0, 1], which is found to be in the form of (ISTT i. As 
seen, (j> in dSTI ) results to be a monotonic increasing function 
of 7. Recalling that 7 takes values in the interval [2, 00), we 
have that the minimum value (pnun is achieved when 7 = 
2. Setting 7 = 2 into the above quartic equation produces 
(/)^ (0^ + 40 — 4) = from which ( |32] | is easily found. 

On the basis of the above results, it follows that Q" (9) 
is zero in Q only for 9 ~ In (p. Moreover, it turns out that 
Q" (9) > for 6* < ln0 and Q" {9) < for ln0 < 6* < 0. 
Then, we may conclude that Q (9) is convex for 9 < In cj) and 
concave for In < 6* < 0. 



The point t/j in ( |35] | belongs to the interval (0, 1] and it is 
such that 

Q(ln V) + Q'(lni/0 (9 - Imp) = Q'(ln V)0 

for Intjj < 9 < 0. This amounts to saying that Q'{lntlj) Intjj = 
(^(In?/)) from which using ( l30l l and ( |63] | we obtain 



The above equation can be easily rewritten as in (|35t in the 
text. 



Appendix E 

The proof of Proposition 2 is divided into the same steps 
used in Appendix C for Proposition 1 . In the sequel, we report 
only the major differences and refer to Appendix C for the 
complete proof. The first step removes the constraints 9n < 
and proceeds computing 9k from 9*j^ is then obtained using 
the following lemma. 



Lemma 3: The solution 9^^ is such that 



=min(0,^K) 



(64) 



where 9x is given by 



9k ~ iT^WK + max ^— . (65) 

e=o,i,...,K-i 



with 



K 



n— 1 n— 1 



(66) 



for j = 1,2, . . . , A and Sq = 0. 

Proof: The proof follows using the same arguments 
adopted in Appendix C and it is not reported for space 
limitations. ■ 
Once 9^ is obtained using ( l64l i and ddSl ), we proceed 
computing 9'^_^. Paralleling the same steps for 9^^, the latter 
is found to be 

9*K_,=min{0,9K-i) (67) 

where 9k-i is obtained solving 

v^^-i -e„ 

s-t. j:Li0n<Sj j = l,2,...,X-l 

J2n=l 0n < Sk - Ok- 

Observing that 5k — 9J^ is always such that 

5k- S*K < 6k-i 
the above problem reduces to the following one: 

mm T,n=i " 

s-t. ELi^«<<5' j = l,...,i^-l 



where 



, ^ , dj for j = 1,2,..., A' -2 

J 1 Sk-9}, ioTj=K-l. 
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The solution 9k-i in (IFTT i has the same form of ( |65] l once K 
is replaced with A' — 1 and the quantities 5i are replaced with 
5'^ for ^ — 1, 2, . . . , A' — 1. Applying repeatedly the same steps 
for any n < AT — 2 leads to the iterative procedure illustrated 
in Algorithm 2 in the text. 
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